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Impurity induced low temperature upturns in both the ab-plane and the c-axis dc-resistivities of 
cuprates in the pseudogap state have been observed in experiments. We provide an explanation of 
this phenomenon by incorporating impurity scattering of the charge carriers within a phenomeno¬ 
logical model proposed by Yang, Rice and Zhang. The scattering between charge carriers and the 
impurity atom is considered within the lowest order Born approximation. Resistivity is calculated 
within Kubo formula using the impurity renormalized spectral functions. Using physical parame¬ 
ters for cuprates, we describe qualitative features of the upturn phenomena and its doping evolution 
that coincides with the experimental hndings. We stress that this effect is largely due to the strong 
electronic correlations. 
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I. INTRODUCTION 

Unlike conventional superconductors following BCS 
mechanism, many aspects of materials showing high tem¬ 
perature superconductivity are still poorly understood. 
Example of such materials include hole doped cuprates 
which show superconductivity at higher temperature. 
These materials show a very rich phase diagram when 
carrier doping and temperature vary. These materials 
not only lack a proper microscopic description of su¬ 
perconducting transition in them, but also their normal 
phases above the transition temperature are poorly un¬ 
derstood. In the over doped regime, electrons in these 
materials behave in a Fermi liquid manner, while a ma¬ 
jor deviation from this behavior is observed at optimal 
doping as well as in the underdoped regime. At low tem¬ 
perature, among many puzzling observations, presence 
of a pseudogap along with the superconducting gap in 
the electronic spectra is observed in this regime^. Micro¬ 
scopic origin of the pseudo-gap is still debatable. How¬ 
ever there is hardly any doubt that many unconventional 
electronic behavior is associated to its presence^. Here we 
focus on understanding some anomalous charge transport 
properties of these materials in the underdoped regime at 
low temperature. 

The effects of pseudo-gap and impurity scattering in 
dc-resistivities of cuprates in the underdoped regime 
is of our concern here. It is well studied in the ex¬ 
perimental literature both in presence and absence of 
impuritie s 345 . Universal linear in temperature behav¬ 
ior of dc-resistivity at optimal doping and its deviation 
both in the underdoped(p ~ log(l/T)) and overdoped 
(p ~ T 2 ) regimes become hallmarks for cuprates. The 
upturn behavior gets enhanced when cuprates are fur¬ 
ther doped with impurities and reminds us a Kondo-like 
mechanism as its origin. Theoretical explanation of it is 
not settled yet and is highly debatable^. However, basic 
elements of theory of such upturn phenomena should con¬ 
tain two elements seriously. The first is the reconstruc¬ 
tion of Fermi-surface and the presence of pseudo-gap in 
electronic density of states due to strong electron-electron 


correlations. This leads to non-Fermi liquid like behav¬ 
ior in the underdoped cuprates. Second is the scattering 
of those non-Fermi liquid quasi-particles with impurity 
atoms at low temperature. We incorporates these two 
important features within a semi-phenomenological the¬ 
ory in the present study. 

Here we limit ourselves to the situations where 
cuprates are doped with non-magnetic impurities only. 
Experimentally, cuprates are further doped with non¬ 
magnetic Zn ions which replaces few Cu ions from Cu- 
O plane and provide a scattering center for charge car¬ 
riers. Doped Zn ions reduces superconducting correla¬ 
tions to certain extent. Also a magnetic field of suitable 
strengths (up-to ~ 50 — 60T) are used to exclude the 
contributions due to superconducting correlations. This 
sets the platform for studying normal state properties 
of cuprates which often shows non-Fermi liquid behav¬ 
ior, particularly in the under-doped regime. It is also 
to be mentioned that, effects similar to the Zn doping 
can be achieved by another method such as electron ir¬ 
radiation with suitable fluence^. In the later case one 
creates point defects in the Cu-0 plane and the gross 
features of the resistivity behavior in this case is almost 
the same as that of Zn substituted one. In such an experi¬ 
mental condition, depending on impurity concentrations, 
an impurity induced supra-linear in T (p ~ T v , v < 
1) to log(l/T) upturn in the resistivity is observed^. 
The later phenomenon hints at localization onset and 
are often attributed to strong electron-electron correla¬ 
tions, spin fluctuation effects, impurity induced electronic 
inhomogeneit y 91011 . Since a detailed microscopic theory 
for underdopped cuprates still lacking, logical direction 
in understanding this material is to use a phenomenologi¬ 
cal theory based on experimental inputs. Recently Yang, 
Rice and Zan g 1312 proposed a phenomenological model 
which can be used to understand properties of electronic 
behavior in the underdoped regime. This theory is quite 
successful in understanding many thermodynamic prop¬ 
erties and ARPES findings^ and also transport proper¬ 
ties like optical conductivity!^ etc. in the recent past. 
We adopt it as a model non-Fermi liquid description suit- 
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able to study low temperature behavior of normal state 
resistivity in cuprates in presence of a pseudo-gap. For 
the rest of our discussion we will call it “YRZ” formal¬ 
ism and in the later section we present a brief review of 
it. In this work, we show that, without adopting addi¬ 
tional inputs one can observe resistivity upturns within 
YRZ formalism. Also the present study reveals the im¬ 
portance of pseudo-gap to this phenomena. 

This paper is organized as follows. In section |TT] we 
first review the YRZ formalism as proposed by Yang, 
Rice and Zan g 1312 . and then derive an expression for 
the dc-resistivity using Kubo formalism and Considering 
the scattering of YRZ-quasiparticles with non-magnetic 
impurities. Then in section cm we numerically evaluate 
the expression for the dc-resistivity and discuss the re¬ 
sults. Finally, in section IIVI we summarize our results 
and conclude. 


II. YRZ-QUASIPARTICLES AND IMPURITY 
SCATTERING 

In the presence of a pseudogap, YRZ formalism is de¬ 
fined with the following ansatz for the coherent part of 
the electronic self energy based on a resonating valence 
bond (RVB) spin liqui d 3312 . 

X yrz Ql,u>,x) = -^ ( 1 ) 

w + e k 

Thus the quasiparticlc propagator or Greens function in 
this formalism takes the form, 


optimal doping and we restrict to the regime x < 0.2. 
In the normal state, the quasiparticle propagator can be 
written in a form, 


£(k,w,x) = ^2 

CK=±1 


9 t(x)Wg(x) 

“-41 1 ) 


( 3 ) 


Above form describes YRZ quasiparticles by the energy 
dispersions, 


K(z) = ^ + «\/(t ) 2 + a ^ (:c) ’ (4) 

and the Luttinger weights are given as, 



Here a = ±1 and = e k ± e£. We consider scattering 
of YRZ quasi-particles with non-magnetic impurities in 
the first Born Approximation. They provide scattering 
centers for quasi-particles and changes their lifetime via 
impurity contribution to the self-energy. Diagrammatic 
representation of the self energy correction due to impu¬ 
rity scattering is shown in Fig[l] Here arrowed parallel 

* ' i'' 

/ = — —^^^— + —s*" 2 ^ ^ + 

dressed YRZ bare YRZ 


G( k,«,*) =- 9t[x) A2 ■ (2) 

W " £k ~ ^ 

The interaction renormalized band-structure disper¬ 
sion is given as = —2t(x)(cos(k x a) + cos (k v a)) — 
At' cos (k x a) cos(k y a) — 2t"(cos(2k x a)+cos(2k y a))— n p (x). 
Here a is the unit cell dimension in the Cu-0 plane. 
Electronic energy band is defined by the renormal¬ 
ized parameters, t{x) = gt{x)to + 3g s (x)Jx/8, t'(x) = 
gt(x)t' 0 and t"(x) = gt(x)t.Q. The Gutzwillcr factors 
which constrain electronic motion in an interacting en¬ 
vironment are given by, g t (x) = ^ and g s (x) = ^+ x yi - 
Here, to, t' 0 , and t'/ are the bare band hopping parame¬ 
ters with t' Q = 0.3toandfg = 0.2f O - J is the magnetic 
energy of the t-J model taken to be J = to/3, x = 0.338 
is the spin susceptibility. e£ = —2t(x) 

(cos(k x a) + cos(k y a)) is a second energy dispersion which 
gives the antiferromagnetic Brillouin zone boundary for 
= 0, which is also referred to as the umklapp sur¬ 
face. The shift in chemical potential g. p {x) is to be 
determined to get the correct hole doping x based on 
the Luttinger sum rule. Self-energy contains a phe¬ 
nomenological input for pseudo-gap energy scale and is 
given as, A pg (x) = A pg {x)(cos(k x a) — cos{k y a))/2, where 
A pg (x) = 3tg(0.2 — x ). In this description x = 0.2 is the 


FIG. 1. Diagram representing renormalized quasiparticle 
propagator due to impurity scattering in the Born approx¬ 
imation. Second and the third diagrams in the right hand 
side represent impurity induced self energy corrections. In 
the present work we consider the lowest order correction, i.e. 
up-to the second diagram in the right hand side. 


lines in the left hand side represents impurity renormal¬ 
ized propagator. The first diagram in the right hand side 
is the bare YRZ propagator as given by Eq[3] and the 
third term represents self energy corrections to the car¬ 
rier propagator due to its interaction with impurities at 
the lowest order. The later is given as^i, 


-‘imp 


(k, iu n ) = Ni ^2 K k - k') | 2 Q( k', iu n ) 


pV 2 -J2S(^',^n) 


= ov 2 i V V 9 t ^ W k' 

V^iUn-EZ 


( 6 ) 


Here impurities potential is assumed to be highly local¬ 
ized i.e. given by v(r) = vS(r). This is justified when 
doped ions sit in the Cu-0 plane and charge carriers 
screen the Coulomb potential. Impurity density Ni/V 
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is given by p and V is the system volume. After per¬ 
forming an analytical continuation we get the frequency 
dependent impurity self energy as, 


(w) = V V-^- 

^ ± V w ~ E ^ +ir i 


-EE 

a:=± k' 




u-EZ 


E E ~ 

a=± k' 


Where Wg = pv 2 g t (x)W£/V. This will be different from 
that of non-interacting case. In case of normal metal, and 
for constant electronic density of states, the real part of 
the Born self energy vanishes, while imaginary part of 
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it becomes constant^. That is not necessarily the case 
here. We see from the above equation that, 

/mSj mp (cu) = 7T EE ^9 E q ) 

cn=± q 

= -pv 2 N YRZ {co)- ( 8 ) 

Here N 1 RZ (co) is the un-normalized density of states for 
YRZ quasiparticles. In previous work s 1516 , where optical 
conductivity of cuprates in the pseudogap state was stud¬ 
ied, 70 = 1 /tq was set to be O.Olfo to regularize the delta- 
function in the expression of the quasiparticle density of 
state s 1417 . Considering impurity scattering within lowest 
order Born approximation, we define impurity renormal¬ 
ized spectral function of quasi-particles as, 


«4(k, co) 


1 


ImY, irnp {u) 


- g t (x) --- 

77 (w - e k - RZ (k : co) - ReT, imp (co)) + {ImY, irnp (u})) 2 


(9) 


The spectral function above includes the incoherent com¬ 
ponent to the bare YRZ spectral function. The incoher¬ 
ent part comes from scattering of the quasi-particles with 
localized impurity atoms. Now we use Kubo formula in 


the linear response regime to numerically calculate the 
conductivity. For a N x N lattice, the expression for con¬ 
ductivity at a finite frequency and temperature is given 
as follows. 


o-(w,T) = 


47re 2 ?i 


(5-5) 

N 2 ca 2 

M-s.-f) 


V,x 


(k) 


dco 


,/V) - /(w + u/) 


CO 


A(k, co')A( k, co + co') 


( 10 ) 


We are concerned about dc-resistivity behavior, we con- at finite temperature, dc conductivity or zero-frequency 

sider the co —> 0 limit of the above expression. Thus conductivity takes the following form, 

I 


CTdc(T) 


4:Tre 2 h 

N 2 ca 2 




do/(-/V)M(kV) 2 - 


( 11 ) 


DC-resistivity p is calculated by taking the inverse of 
&dc{T). Here the frequency derivative of the Fermi-Dirac 
distribution function is given as, 

/,( “>=L (^ihd=iL{l - T snh (£)) 

=-j? sech! (£) ■ (i2) 

This function peaks near co = 0 and has width 2 T(kg = 
1). Thus for numerical calculations, we can fix our co- 
cutoff in calculating dc-conductivity from Eallll as few 


T. 

C-axis conductivity. In Cuprates, charge carriers inter¬ 
act in the Cu-0 or ab-plane and studies related to these 
materials are limited to two dimensional model. However 
in a real material, such Cu-0 planes stack along c-axis 
and carriers can tunnel between two neighboring planes 
through tunneling mechanism. Conduction along c-axis 
is also studied extensively^. In the case of non-magnetic 
impurity doped samples, an upturn in the low tempera¬ 
ture regime is also observed^. We revisit the same here 
within our formalism. To calculate conductivity along c- 















4 


axis, we follow Chakravarty et al^ and introduce tunnel¬ 
ing matrix element between adjacent Cu-0 layers which 
will replace v x in the earlier expressions for conductivity 
in EafTTl with, 


f_i_(k)c = (cos (ak x ) — cos (ak y )) 2 . (13) 


Here t± is the tunneling matrix element for quasiparti¬ 
cles in cuprates along the c-axis. It is usually in the range 
0.1 — 0.15eV. The parameters c and a are the lattice di¬ 
mensions along the c axis and a or b-axis. With these 
descriptions, we conclude our formalism. 


III. RESULTS 

In this section we discuss the findings of our formal¬ 
ism. We show plots of various quantities for 100 x 100 
lattice and for to = 1 eV. Other parameters are taken 
from previous works by one of the present author s 1516 . 
Coefficient of Eq[ll]is set such that p a b{T ) is given in the 
/if2 .cm unit. Before presenting the results, we estimate 
the validity of our approximation. We consider the lowest 
order correction to the quasiparticle self energy that gives 
a frequency dependent contribution and limit ourselves 
to the 2nd diagram in the right hand side of the Fig jT] 
Next order correction, as shown in the diagram next to 
it, is given as pv 3 (IV 1 RZ ^j while the previous diagram 
gives a pv 2 ( N YRZ ) contribution^!. Thus our perturba¬ 
tive approach depends on the condition that the ratio of 
this two successive contributions, i.e. vxN YRZ « 1. In 
our system we see that near the Fermi level, N 3 RZ ~ 0.1 
per eV and if we assume the energy scale of the coupling 
between impurities and quasiparticles v ~ 100 meV, then 
v x N yrz ~ 0.01 << 1 and the condition is well satisfied. 
Within this approximation, impurity induced correction 
to the self energy is determined by a parameter pv 2 than 
v alone. We consider it as a free parameter and vary it 
within certain range assuming that earlier conditions are 
satisfied. 

In Fig[5] we show plots of momentum integrated un¬ 
normalized density of states at x = 0.05 as a function 
of energy for various impurity strengths. It is observed 
that density of states has a valley near u> ~ p and two 
peaks are separated by twice the pseudo-gap energy scale, 
i.e. 2A° However since the pseudogap is not constant 
for all momentum states and contains nodes and anti¬ 
nodes in it’s structure, it affects the density of states 
profile within the range (—A 6 S , A° ). Under its influence 
density of states near the Fermi surface never vanishes, 
rather becomes energy dependent. This is one aspect 
where the pseudogap physics becomes different from that 
of the simple metals. The plot for pv 2 = 0.00 corresponds 
to the bare quasiparticle case where a regularization in 
the Dirac delta function was made by introducing an ar¬ 
tificial width 70 = O.Olto- In this figure we also plot the 
density of states for various impurity strengths and con¬ 
firm that coherence in the quasiparticle density of states 


survives within our choices of impurity strengths. In the 
inset of the same figure we show the behavior of the den¬ 
sity of states near the Fermi surface within the energy 
scale of our interest. It clear that the unlike the case of 
metal, density of states is strongly energy dependent in 
this regime as expected in the pseudogap regime. 



FIG. 2. Un-normalised density of states of quasi-particles vs. 
frequency in eV at different impurity strengths. Here hole 
concentration x=0.05. A non-constant nature of the density 
of states near the Fermi surface is shown in the inset. 

From our experiences in calculating transport coeffi¬ 
cients in metals, it is tempting to relate the density of 
states to the scattering rate of the quasipartic.les or the 
resistivity. However from Eq[lT] and EqJT 2 ] we see that 
it is not exactly the electronic density of states near the 
Fermi level(w = 0) within the energy scale ~ T (~ few 
meV) determines the scattering rate of the quasiparticles 
and hence the resistivity. A quantity which we can name 
transport density of states or tDOS= )Ck ^A(k, w) 2 , as 
opposed to the density of states obtained from the spec¬ 
troscopic experiments and is given by ]U) k A(k, oj) de¬ 
termines the scattering rates and hence resistivity. In 
Fig|3] we see that the density of states near the Fermi 
level within an energy interval of few meV is strongly en¬ 
ergy dependent. Also, with increasing impurity strength, 
the peaks in the density of states profile decrease and 
the spectral weights are transfered towards the Fermi 
level. At stronger impurity strength, the curvature of 
the density of states profile near oj = 0 changes. Since 
A(k, oj) is positive for all k, all these features, though 
in a qualitatively different way, should manifest in the 
tDOS also. Now for qualitative understanding, we can 
assume tDOS ~ (a — 6 sgn(w)|a;| a ), where a, b > 0 and 
a are constants. Using Eq[lTl we see a ~ f^ A du>(a — 

6 sgn(w)|o;| Q ) sech 2 (w/T). The thermal factor is positive 
for all the values of oj and behaves like a Gaussian whose 
width is given by T. Now scaling oj with temperature, we 
get I = cr ~ dx(aT — 6 T 1 +a sgn(x)|x|“) sech 2 (x) = 
AT — bt 1+ “ in the limit T « A. Here a is deter¬ 
mines the curvature of the tDOS profile and along with 
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B, determines the deviation from the uniformity. A more 
appropriate fitting of tDOS should contain other powers 
of u> and they will also contribute to the resistivity. It 
is also nice to comment here that, if for some interac¬ 
tions tDOS of quasi-particles fits with a = —2, resistiv¬ 
ity becomes linear in T. Here all the constants have de¬ 
pendence on the pseudogap, hence doping as well as the 
impurity scattering. In our numerical calculations we see 
that the competition between the such terms as observed 
in the above expression, give rise to the upturn. The de¬ 
pendence on the velocity v x in tDOS is also reflected in 
qualitative differences between resistivity profiles in ab- 
plane and along c-axis, as is shown in our results. We 



FIG. 3. Resistivity in the ab-plane vs. temperature at differ¬ 
ent impurity strengths at a fixed doping x=0.05. An upturn is 
seen at pv 2 = 0.15 around 10K which gets stronger in higher 
impurity strengths. 



FIG. 4. Inter plane or c-axis resistivity vs temperature at 
different impurity strengths at a fixed doping x=0.05. Here 
also we see upturn behavior stronger than that of the ab-plane 
and also gets stronger with the increase in impurity strengths. 

present the plots for resistivity with temperature at dif¬ 
ferent impurity strengths and doping and focus on the 
effects of pseudo-gap phase which is most prominent at 
low doping. We first consider x = 0.05, and look for 



FIG. 5. The ab-plane resistivity vs. temperature with T— axis 
in a logarithmic scale at different impurity strengths. Hole 
doping is fixed at x=0.05. 

resistivity vs. temperature behavior at different impu¬ 
rity strengths. In Figj3] and FigU we see weak upturn 
for both p a b and p c at low temperature, which increases 
with increasing impurity strength. To study the appear¬ 
ance of a log(l/T) behavior in resistivity associated with 
this upturn, we re-plot the Fig. [3] with temperature axis 
in a logarithmic scale and is shown in Figj5] We see that 
at higher impurity strength, resistivity profile in the up¬ 
turn regime moves toward a log(l/T) one. This is in 
accord with experimental findings^. Next, to study 



FIG. 6. Resistivity in the ab-plane vs. temperature at a fixed 
value of pv 2 = 0.15 at different hole concentrations. Resistiv¬ 
ity upturn is seen to become weaker with the increase in hole 
concentrations i.e. as one moves away from the pseudogap 
phase. 

the effects of carrier concentration or pseudo-gap in the 
upturn formation we plot again resistivity vs tempera¬ 
ture at different doping with fixed impurity strength. We 
choose pv 2 = 0.15, as an upturn is visible at this impu¬ 
rity strength. The ab-plane and the c-axis resistivities 
are plotted for various doping in FigjG] and Fig. [7] re¬ 
spectively. We see that, in both cases, the upturns get 
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FIG. 7. Inter plane or c-axis resistivity vs temperature at a 
fixed value of pv 2 = 0.15 at different hole concentrations. It 
also shows weaker upturn with the increase in hole concen¬ 
trations as found in the case of ab-plane resistivity. 

weaker with increasing hole concentrations, i.e. as one 
moves away from the pseudo-gap phase. 

IV. SUMMARY AND CONCLUSIONS 

The main observations from this work are as follows. 
Unlike the case of metals, impurity contribution to the 
quasi-particle self energy in the pseudo-gap phase of 
cuprates is frequency dependent and is determined by the 
density of states profile. Replacing the constant width 
of the single particle states by this frequency dependent 
self-energy, can lead to an upturn in the temperature 
variation of the resistivity at low temperature which is 
often considered as precursor to a metal-insulator tran¬ 
sition. Also with higher impurity strengths, i.e. higher 
values of pv , we see a proximity towards log(l/T) like 
behavior in resistivity. Variation of the upturn with car¬ 
rier concentration is qualitatively in accord with experi¬ 
ments. Now we compare our formalism with previously 


proposed theories. We observe, emphasis on either or 
a combination of the 1) Kondo like scenarios, 2) weak 
localization corrections, 3) impurity induced inhomoge¬ 
neous electronic distribution picture and 4) pseudogap 
and strong electron-electron correlations in all the pro¬ 
posed theories depending on the regime of their interest 
in the cuprate phase diagra m 242526 . For example authors 
in Refi 2526 focus on optimally and slightly underdoped 
cuprates and develop their theory and studies the disor¬ 
der contribution to the upturn phenomenon in a metalic 
Fermi liquid background. In their studies they limit elec¬ 
tronic correlations within weak interaction regime while 
disorder contribution is calculated within a numerically 
exact way. They find importance of impurity induced 
magnetic droplets in the formation of resistivity upturn. 
Though such a picture may be valid near optimal dop¬ 
ing, but certainly not in the deep underdoped regime 
where strong electron-electron correlations is of utmost 
interest. We focus on the later regime. Since there is no 
well accepted microscopic description valid in this regime, 
we adopt a simpler phenomenological picture, namely 
YRZ-theory to incorporate non-perturbative aspects of 
electron-electron correlations. Thus we show that even if 
we limit ourselves to the perturbative Born approxima¬ 
tion to consider the elastic scattering between strongly 
correlated electrons, the upturn behavior in resistivity 
can appear at low temperature. This is certainly a no¬ 
table observation which puts more weight to the elec¬ 
tronic correlations than the disorder effects. From our 
analysis, we conclude that the presence of the pseudo¬ 
gap is enhancing the upturn phenomena which is also ob¬ 
served in the previous theoretical proposal^. Observed 
log(l/T) like behavior of resistivity at higher impurity 
concentrations is also a prominent experimental feature 
in these materials and its origin is highly debatable^. 
However, it needs to be mentioned that in this work, we 
limit ourselves to isotropic, zero range, non-magnetic im¬ 
purities. More realistic impurities providing momentum- 
dependent impurity potential^ and impurities with net 
magnetic moment, e.g. Ni doping^ are surely worth 
studying. These are left for future studies. 
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